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ABSTRACT 

Context. Dust is a good tracer of cold dark clouds but its column density is difficult to quantify. 

Aims. We want to check whether the far-infrared and submillimeter high-resolution data from Herschel SPIRE and PACS cameras 
combined with ground-based telescope bolometers allow us to retrieve the whole dust content of cold dark clouds. 

Methods. We compare far-infrared and submillimeter emission across L183 to the 8 pm absorption map from Spitzer data and fit 
modified blackbody functions towards three different positions. 

Results. We find that none of the Herschel SPIRE channels follow the cold dust profile seen in absorption. Even the ground-based 
submillimeter telescope observations, although more closely following the absorption profile, cannot help to characterize the cold 
dust without external information such as the dust column density itself. The difference in dust opacity can reach up to a factor of 3 
in prestellar cores of high extinction. 

Conclusions. In dark clouds, the amount of very cold dust cannot be measured from its emission alone. In particular, studies of dark 
clouds based only on Herschel data can miss a large fraction of the dust content. This has an impact on core and filament density 
profiles, masse and stability estimates. 
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1. Introduction 

Dark clouds are the places where stars form. We try to follow 
the different steps that lead from a low-density cloud to main- 
sequence star. Some steps have been clearly identified like the 
Class 0 to Class III protostar evolution fLada|1987 Andre et al. 

11993) but, the early phases, such as the (formation of a prestel¬ 
lar core (PSC), the collapse of this core, are still not well un¬ 
derstood. Studies of clouds and of PSCs first attempt to deter¬ 
mine their mass. This is not a simple task. The first difficulty 
is to know the distance of the object, and the second is to trace 
the column density of the material itself. The main components, 
H2 and He are not directly visible, with the exception of a lit¬ 
tle Ht seen in absorption in the UV at the cloud edges, and 
surrogate tracers are needed. Molecules are not good tracers in 
general because CO, which is the standard tracer, is depleted in 
the P SCs (|Lemme et al.|1995| Willacy et al. 11998) Tafalla et al. 


|2002| Pagani et al. |2005[|Brady Ford & Shirley||26ll| ). Tracers 

of the PSCs like NH3 or NtFF do not extend beyond the PSCs 
themselves, and contrary to CO, their peak abundance is vari¬ 
able. Detailed radiative transfer models are therefore needed to 
retrieve the H2 + He densities via the modelling of their colli¬ 
sions with the tracers (followed by integration along the line of 
sight to obtain the total column density). On the other hand, dust 
is traceable from the edge of the clouds to the centre. Its relative 
abundance to H2 is not accurately known but is thought to be in 
the range ~1/130—1/100 in the Milky Way (Flower et al. 2005 
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Compieg ne et al.|2011). It c an be traced either by its extinction 
of background stars ( |Wolf||1923| Bok 1956| Bok & Cord well 
1973) with a higher efficiency in the near-infrared (NIR, Lada 
etal.|1994| |Lombardi & Alves|200 1 [ |Lombardi|2009 1 or of back¬ 


ground diffuse light in the mid-infrared (MIR, Bacmann et al. 


the b 


2000), by its scattering of the interstellar radiation field (ISRF) in 
IR and MIR (Lehtinen & Matti la| 1 996[|Juvela et al.|2006[ 


Lefevre et al.|2014 1, or by its emission in the far-infrared (FIR) 

(Ward-Thompson et al.: 1 994) . However, all these methods have 
difficulties. The use of background stars do not allow to reach a 
high spatial resolution outside the galactic plane and become ab¬ 
sent typically when Ay > 40 mag. Scattering is a promising but 
difficult method, highly dependent on the type of grains (size dis¬ 
tribution, extinction, albedo, phase function), on the background 
intensity, and on the ISRF anisotropy and strength (see Lefevre 


et al. 


2014 for an exhaustive study). Emission depends on the 


know edge of the grain properties, spectral index and tempera¬ 
ture (e.g. Juvela & Ysard 2012| for a discussion of the degeneracy 
of the problem). The advantage of dust emission measurements 
in the FIR and submillimetre domains is that it remains optically 
thin up to N(H 2 ) ~2 x I0 24 cm 2 at 300 pm and even higher at 
longer wavelengths. Dust observed in emission is therefore able 
to trace the whole cloud content from the edge to the centre. 
However, it remains difficult to convert this emission into the 
actual column density of dust. 

Before the launch of Herschel ( {Pilbratt et al.|2010| ), FIR mea¬ 
surements of dust emission in the range 200-500 pm were scarce 
(e.g. Stepnik et al.|2003 1 and many studies limited themselves to 
fit the dust in emission by building single spectral energy dis- 
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Fig. 1 . LI83 Herschel/SPIRE map at 250 fim. Contour levels ev¬ 
ery 20 MJy sr 1 . The cut through the region is traced as a ~6° 
tilted line with respect to the vertical (equatorial North) direc¬ 
tion. 



tribution (SED) fits based on a mix of measurements shorter 
than 200 gm and longer than 850 fim. However, Pagani et al. 


( j2004| l showed that in a dark cloud without embedded sources, 
the 200 fim emission traces only the envelope of the cloud (A v 
< 7.5 mag from the surface) and totally misses the bulk of the 
dust too cold to emit significantly at that wavelength. Hence, it 
could be deduced from that result that single SED fits were not 
realistic because short wavelengths are mostly tracing emission 
from the envelope and long wavelengths emission from the core. 
There was no connection between the two sides of the (modi¬ 
fied) blackbody curve and therefore no physical meaning to that 
fit. The question arises again today now that all wavelengths 
between 100 and 1200 yum have been sampled in a number of 
clouds with several telescopes, including Spitzer ([Werner et al. 


|2004j l, AKARI ( [Murakami et al.| 2007},gr ound-based bolome¬ 
ters, and above all, Herschel dPilbratt et al.||2010[ > and Planck 
(Taub er et al.|2010| , to find out if we can accurately retrieve 
the dust content from fitting SEDs alone. Indeed, dust at 10 K 
is brightest at 300 yum, while dust as cold as 6 K is brightest 
at 500 fim. Both wavelengths are inside the range observed by 
SPIRE, the FIR camera of Herschel ( [Griffin et al.||2010| and 
not far from the two highest frequency channels of PlanckfHFl 


Lamarre et al. 

2010). Using Herschel data, Nielbock et al. 

2012 ) achiever 

a successful temperature analysis of B68 but 


they combined data from NIR to submm. Roy et al. (2014 1 at¬ 
tempted to fit two cores (B68 and L1689B) using Herschel data 
alone. They both find temperatures lower than provided by single 
blackbody fittings of the SEDs. Some aspects of dust emission 
fitting of mock clouds have also been numerically explored by 
|Malinen et al.| ( [20TT[ ). 


2. Observations 


This work is based on maps of L183 taken at 8, 100, 160, 250, 
350, 500, 870, and 1200 fim. The 8 yum data were taken with 
Spitzer/IRAC, [Fazio et ~aT7||2004[ published in |Steinacker et al. 
201 Oj i. Herschel/PACS took the 100 and 160 yum maps, and 


Herschel /SPIRE the 250, 350, and 500 /jm maps (Pilbratt et al. 
|2010||Poglitsch et al.|2010[|Griffin et al.|2010[ ). These maps are 
partly published in Juvela et al. ( 2012 for the 250 yum data) and 
will be published in Montillaud et al. (2014 >; Juvela et aI7| ( 2014] >, 
and Lefevre et al. (in prep.). Lefevre et al. will also publish two 
new APEX/LABOCA and IRAM-30m/MAMBO maps at 870 
and 1200 fin l, respectively. In this Letter, we only present a cut 
inclined eastward by ~6° with respect to the vertical (Fig. |T|). It 


is adjusted to go through the main prestellar core (PSC) and the 
northern prestellar core (N-PSC, Figs. [I]&[2]). The minimum in 
extinction between the two cores is named the saddle. Details 
about data reduction will be exposed in Lefevre et al. (in prep.). 
Shortly, Herschel data acquired as part of the Galactic cold 
cores key project (Juvela et al.|2010l were re-reduced with the 
Herschel interactive processing environment (HIPE) v. 12.1, us¬ 
ing the official pipelin e ([Ott|2010 1 , including Scanamorphos v.23 
for PACS data (Roussel [201 3|). They were subsequently colour- 


corrected, zeroed, and rescaled using IRIS (Miville-Deschenes 
& Lagache 2005) and Planck data (Juvela-subm). The calibra¬ 


tion error is a small fraction of the background subtraction uncer¬ 
tainty which we estimate at 10% for SPIRE and 15% for PACS. 
We reduced LABOCA data using BoAQ and MAMBO data us¬ 
ing MOPSIC0 Both LABOCA and MAMBO data lack the large 
scale, low surface brightness part of the cloud because of the 
observing method of MAMBO (in-source OFF subtraction) and 
the data reduction method for LABOCA, which requires strong 
filtering to remove sky fluctuations. To recover most of the lost 
signal, we combined these data with Planck data that are colour- 
corrected and interpolated to match the LABOCA and MAMBO 
filter band-passes. For MAMBO, the correction is no more than 
10%, while for LABOCA the correction amounts to 30 %. The 
final uncertainty is estimated to be ~20% for both. 


3. Analysis 

Figure [2] shows the extinction profile along the cut (as defined 
in Fig.jTj), derived from the 8 fim data degraded to the resolu¬ 
tion of, and compared with, the 250 fim and the 870 fim cuts. 
The 8 fim opacity values are only a qualitative representation 
of the expected extinction. It is based on a previous estimate by 
Pagani et ak| (|2004 1 , on the necessity to reach zero extinction 


at the edges of the cloud, and on the estimated ratio in column 
density between the PSC and the N-PSC, which is 1.5-2 at the 
MAMBO 12" resolution (after smoothing to 18", the ratio is in 
the range 1.4-1.7). The extinction is represented as a range be¬ 
tween these two ratio values. These values are only indicative 
since the 8 yum absorption map is in fact strongly contaminated 
by scattered light and cannot be safely converted to an opacity 
map without a 3D model taking the scattering due to micron¬ 
sized grains into account ( |Steinacker et al.|[2010| Pagani et al. 
|2010[ [Lefevre et al.||2014) . A better estimate of the extinction 
based on NIR and MIR absorption and scattering, including the 
8 fim map, will be presented in Lefevre et al. (in prep.). 

While the strong features from the 8 yum data and the 870 fim 
data are clearly correlated, the weak peaks southward (negative 
offsets in Fig. [2J seen at 8 fim are missing at 870 fim. This is due 
to the loss of extended features with LABOCA. These features 
are too small to be resolved by Planck and are detected only as 
a smooth increase in emission in this Planck +LABOCA com¬ 
bined data cut. The 250 fim emission detects the two southern 
peaks (with a small shift for the most opaque one, most proba¬ 
bly due to an anisotropic heating). However, it does not correctly 
trace the main PSC since the main peak and the northern peak 
have almost the same intensity instead of the ratio 1.4-1.7 we 
expect from submm observations. Both peaks are also too weak 
compared to the saddle intensity. This is reminiscent of the sim¬ 


ilar result reported by Pagani et al. (20041 in which they showed 


that the ISOPHOT 200 yum map does not follow the dust col¬ 
umn density in the inner part of the cloud. The low resolution 


1 http://www.eso.org/sci/activities/apexsv/labocasv.html 

2 http://www.iram.es/IRAMES/mainWiki/CookbookMopsic 


2 





















































































L. Pagani et al.: Can we trace very cold dust from its emission alone ? 




Offset (arcminutes) 


Fig. 2. L183 cut. The 8 gm data have been smoothed to 18" 
resolution to match both the Herschel/SPIKE 250 fi m and 
APEX/LABOCA 870 gm common resolution. The 8 fim absorp¬ 
tion cut is converted to relative extinction (see text). The grey 
filling represents the possible range of 8 gm ratio between main 
and northern peaks. All curves are aligned at the saddle point (at 
+2'33"offset). 

(90") of the ISOPHOT map could have been a partial explana¬ 
tion, but it can now be seen that resolution is not an issue and the 
dust is simply too cold to contribute significantly to the 250 fim 
emission in the central parts of the cloud. Indeed, the peak inten¬ 
sity ratio between the southern peak at offset -4'40" and the two 
main peaks is about 2 at 250 fin t, while it is larger than 3 for the 
northern peak and in the range 4-5 for the main PSC in terms 
of relative opacity at 8 fin t. Figure [3] shows the dust emission 
along the cut at wavelengths from 100 to 1200 fim. The resolu¬ 
tion of all the data has been aligned on that of the SPIRE 500 fim 
channel (37"). The longer the wavelength, the better the tracing 
of the two dust peaks. Figure [4] displays the SEDs for the three 
points of interest along the cut: the main PSC, the saddle, and 
the northern PSC. For all three positions, we tried to fit the SED 
with either two or three modified blackbodies, one at 17K, one 
at ~ 10 K, and an optional one at 6 K, 

n n 

I v = y,B v (T d ,,)r Vo (-f = V B v (T d , i > VoAi m H N(H 2 , i )(-f ,(1) 
vo v 0 


i=l 


Fig. 3. L183 cut for dust emission. This cut includes data 
from Herschel PACS and SPIRE, APEX/LABOCA, and IRAM- 
30m/MAMBO data, all convolved to 37" resolution, shown one 
by one. The fluxes are scaled down by the amount indicated to 
the right of the wavelengths to align all the fluxes on the saddle 
point (except at 100 fim). The 8 fim opacity range is displayed 
(black lines filled in grey), the 250 fim is traced in blue for easier 
comparison with Fig.pl Click here to change wavelength . The 
FIR wavelengths (except the 100 fim one) are shown together in 
the Appendix (Fig. |A.1[) 


with n = 2 or 3. The parameter B v is the Planck function at dust 
temperature T d , vo = 1 THz (A = 300 fim), fi = 2.33 the mean 
molecular weight, mn = 1.67 x 10 ~ 24 g the proton weight, N(H 2 ) 
the gas column density, and k Yo the dust opacity at 300 fim. The 
= 0.111 cm 2 g _1 


opacity at 300 fim varies from 0.02 +0.07 (n = 2) to 0.058 (n = 
3), almost a factor of three, in the case presented in Fig. [4} and 
could go up to 0.1 (a factor of five higher) in the extreme case 
where /i f)K ~ 3.9. As expected the difference is less for the sad¬ 
dle and the northern PSC but the increase in opacity can still be 
> 50%. There is no way to discriminate between the fits simply 
based on these data. Of course, three blackbody fitting instead of 
two is still far from the real case (a continuous variation of tem¬ 
perature from ~ 17 K at the surface of the clouds to ~ 6 K in the 
densest cores, [Zucconi et al.|2001[ [Evans et al.|2001| > but comes 
much closer to it already. 

The fits we find with an imposed 6 K blackbody are not 
unique. Since a solution can be found without a 6 K compo¬ 
nent, the opacity of this component can be varied from 0 up to 
0.1 depending on the other parameters, T dust and ft in particular, 
with similar ^ 2 values. Constraints must come from other obser¬ 
vations, which we discuss now. 

For the main PSC, N(H 2 ) is in the range 4.5+1.7 x 10 22 


i+i.o. 
- 0.8 ' 


V23 


at density 1 x 10 6 cm 3 from 

Ossenkopf & Henning 

(1994 

) for a 

gas-to-dust ratio of 133 ( 

Compiegne et al. 2011), which is com- 

patible with Ysard et al. 

(2013). The spectral index f3 describes 


encompasses the column density determined via gas modelling. 
Pagani etal. (2007) indicate N(H 2 ) ~1.0x 10 23 cm 2 from N 2 H + 


the modification of the dust opacity k v with frequency. The main 
blackbody is optimized for all three parameters, Tj,/), N(H 2 ). 
The cold blackbody, if included, is set at 6 K, and f3 and N(H 2 ) 
can be adjusted. The third blackbody at 17 K, with [3 — 1.8, is set 
to fit the grains on the cloud surface ,which contribute at 100 fim. 
Their only free parameter is N(H 2 ). Their temperature and spec¬ 
tral index are typical of diffuse and cloud surface dust temper¬ 
atures ( |Zucconi et al.|[200T| jPlanck Collaboration|[2011| l. Their 
contribution is always less than 0.3% of the total mass and is not 
discussed any further. We independently adjust fi values in the 
range 1.5-4 with the constraint that /f T | < /j K - The fits are op¬ 
timized so that x 1 < 0.5 in all cases (all parameter values are in 
Fig . [4]). The striking result is that the fit is just as good with n = 2 
(blackbodies at ~ 10 and 17 K) or n = 3 ( 6 , ~ 10, and 17 K). The 
X 1 values remain basically identical. Towards the main PSC, the 


data when averaged in a 37" beam (after correcting for the new 
N 2 H + -H 2 collisional coefficients, Lique et al. 2015). However 
N 2 H + only traces the dense region and not the envelope of the 
cloud. From C ls O measurements (Pagani et al. 2005), we infer 
a cloud envelope column density of ~2 x 10“ cm~ 2 . The total is 
~1.2 x 10 23 cm -2 . The difference with the n = 2 fit (ratio of 2.5) 
is larger than the different uncertainties involved here. With the 
new N 2 H + -H 2 collisional coefficients, we also found that the gas 
temperature in the core is 6 K (7 K in [Pagani et al. 12007} where 
the density is > 5 x 10 5 cm~ 3 , which is a density high enough to 
efficiently thermalise gas and dust (n > lx 10 5 cm ~ 3 is required, 
|Goldsmith[2001) . Therefore, dust and gas temperatures and col¬ 
umn densities can be made consistent only if we introduce a 6 K 
dust component in the fit. Based only on Herschel data,[Roy et al. 
( 2014) also find a temperature lower than that given by a single 
SED fit (9.8 K instead of 11.6 K in the case of L1689B at its cen¬ 
tre), and their fit is more realistic than the one we present here by 
using a continuously varying temperature with radius. However, 
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Fig. 4. SEDs of the three points of interest defined in Fig. 
[2] along the cut (in apparition order, the main PSC, the 
saddle and the northern PSC). Top row, the SEDs are 
fitted with three modified blackbodies, bottom row, with 
two modified blackbodies. Opacity is given at 300 jjm. 
[Click to display the three cases in turn . The three plots are 
shown separately in the Appendix (Fig.|A.2|. 


despite this more sophisticated analysis, Roy et al.’s method fails 
to reveal the very cold dust (6-8 K) in the heart of L1689B that 
we have identified from NiH + measurements (Bacmann et al., in 
prep). These measurements indicate a 37" beam-averaged col¬ 
umn density of 7-11 X 10 22 cm -2 , again a factor 2 to 3 times 
higher than Roy et al.’s 3.5 x 10 22 cm -2 , which does not even 
take the contribution from the envelope devoid of N 2 I-E emis¬ 
sion into account. Another attempt by Marsh et al. (2014) gives 
better results. They trace cold dust down to 6 K but they use a 
priori density and temperature profiles and therefore their results 
somewhat depend on their input parameters in contrast with Roy 
et al. approach. 

It is clear from these two cases that very cold dust exists and 
is not identified as such solely by its emission. This is because 
of the well-known fact that warmer dust (even as low as 10 K) 
outshines the very cold dust because of the non-linearity of the 
blackbody function (at 6 K, the Rayleigh-Jeans approximation is 
not valid above 100 GHz). This is clearly visible in Fig. [4] where 
the 6K blackbody contribution is smaller than the ~ 10 K con¬ 
tribution, even at A = 1.2 mm and even when it represents 70% 
of the total mass (45% of the signal in that case). Though FI83 
and F1689B could be thought of as peculiar cases, we believe 
they are only clear illustrations of a general problem, since dust 
below 9 K is predicted to occur for radial extinctions as low as 
Ay= 10 mag (Zucconi et a l.|2001 1 , in standard ISRF conditions 
(no local heating source). It is clear that studies of filaments, and 


cold dark clouds in general, based solely on dust emission (and 
worse, single component SED fits), even including Herschel or 
Planck data, should meet the degeneracy problem we expose 
here and will miss a large fraction of the mass in these objects. 


Even more sophisticated modelling, as presented by Roy et al. 
( 2014} , misses a large fraction of the mass. In the present two 
cases, F183 and F1689B, 30 to 70% of the total mass is lack¬ 
ing. The dust column density must be retrieved by other means 
(dust in extinction, molecular emission). Therefore, critical and 
threshold masses, core and filament stabilities, and density pro¬ 
files in studies limited to dust emission should be considered 
with caution. 
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Appendix A: Multi-layered picture 

Figure |AT| shows the layers of Fig.[3]assembled together for di¬ 
rect comparison, for pdf viewers not understanding Javascript, 
and for printers. Similarly, Fig. |A.2| shows the layers of Fig. [4] 
(the three positions tested with two or three blackbody compo¬ 
nents) side by side. 
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Offset (arcminutes) 

Fig. A.l. LI83 cut for dust emission. It includes data from Herschel PACS and SPIRE, APEX/LABOCA, and IRAM-30m/MAMBO 
data, all convolved to 37" resolution. The fluxes are scaled down by the amount indicated to the right of the wavelengths to align 
all the fluxes on the saddle point. The 8 pm opacity range is displayed (dashed lines filled in grey), the 100 pm cut is omitted since 
it only traces the cloud content at its surface. This is the developed version of Fig. [3] The colour is changed for each wavelength to 
help separate them. 
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Fig.A.2. SEDs of the three points of interest defined in Fig. [ 2 ] along the cut (from left to right, the main PSC, the saddle and 
the northern PSC). Top row, the SEDs are fitted with three modified blackbodies, bottom row, they are fitted with two modified 
blackbodies. Opacity is given at 300 pm. This is the developed version of Fig. [4] 


















